Energy as an Entanglement Witness for Quantum Many-Body Systems 



Mark R. Dowling0 Andrew C. Doherty, and Stephen D. Bartlett 
School of Physical Sciences, The University of Queensland, Queensland 4072, Australia 

(Dated: 4 January 2005) 

We investigate quantum many-body systems where all low-energy states are entangled. As a 
tool for quantifying such systems, we introduce the concept of the entanglement gap, which is the 
difference in energy between the ground-state energy and the minimum energy that a separable (un- 
entangled) state may attain. If the energy of the system lies within the entanglement gap, the state 
of the system is guaranteed to be entangled. We find Hamiltonians that have the largest possible 
entanglement gap; for a system consisting of two interacting spin-1/2 subsystems, the Heisenberg 
antiferromagnet is one such example. We also introduce a related concept, the entanglement-gap 
temperature: the temperature below which the thermal state is certainly entangled, as witnessed by 
its energy. We give an example of a bipartite Hamiltonian with an arbitrarily high entanglement-gap 
temperature for fixed total energy range. For bipartite spin lattices we prove a theorem demon- 
strating that the entanglement gap necessarily decreases as the coordination number is increased. 
We investigate frustrated lattices and quantum phase transitions as physical phenomena that affect 
the entanglement gap. 
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I. INTRODUCTION 

Understanding and quantifying the properties of quan- 
tum many-body systems is a central goal of theoretical 
condensed matter physics. Progress is often hindered by 
an incomplete understanding of the highly non-classical 
entangled states that occur naturally as the ground and 
thermal states of many systems. Entanglement is per- 
haps the most counter-intuitive feature of quantum me- 
chanics, and results in stronger correlations than can be 
present in any classical system 0, . Recently, entangle- 
ment has been recognized as an important resource in the 
emerging field of quantum information science 3] , which 
has led to new tools that may enhance our understand- 
ing of the role of entanglement in quantum many-body 
systems. 

Much recent work has focused on quantifying the en- 
tanglement naturally present in the ground state of stan- 
dard models of coupled quantum systems, particularly 
spin chains. In 0, H, IE S IS IS the role of entangle- 
ment in a quantum phase transition [TH is investigated. 
In one-dimensional chains, the amount of entanglement 
between a length of spins and the rest of the chain ap- 
pears to depend only on the universality class of the 
model at the phase transition |S 0- Various quanti- 
ties associated with entanglement have been shown to 
display universal scaling behavior at phase transitions in 
one dimension 0, |S U- Also, it appears that proper- 
ties of entanglement between spins, such as the entangle- 
ment length defined in [l^ . are sometimes able to char- 
acterize phases of the system better than any correlation 
length |l|. 

Restricting to many-body systems where each system 
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interacts with only a finite local neighborhood (which 
we refer to as local interactions) very strongly constrains 
the quantum states that must be considered. For ex- 
ample, there exist quantum states that are far from the 
ground state of any local- interaction Hamiltonian 
For finite systems with local interactions and an energy 
gap it was shown in [l^ that there are strong bounds 
on both the correlations and the entanglement in the 
ground state. The fact that local interactions strongly 
limit the entanglement that can occur for Hamiltonian 
systems with local interactions on a line or a plane has 
been used to develop new approximation schemes for sim- 
ulating quantum dynamics limilllliSiSill. There 
is now a large literature on the entanglement properties 
of the ground states of Hamiltonian systems; we refer the 
reader to [HililiSlllIll and references therein. 

Although the ground state plays an essential role in 
understanding physical systems, at finite temperature it 
is the thermal state that is of the greatest interest. The 
nature of entanglement in the thermal state of condensed 
matter systems was first studied by Nielsen |2S| , who in- 
vestigated how entanglement in the thermal state varied 
with temperature and other parameters of simple systems 
consisting of two coupled spins. Subsequent work has 
investigated similar questions for quantum many-body 
systems HIMIIllESIlSilllil. A recent exper- 
iment demonstrates that entanglement can affect ther- 
modynamic properties of a system at high temperature 

Thus, it seems that many physical phenomena involv- 
ing just the ground or thermal states in condensed matter 
systems may be associated with the nature of entangle- 
ment in the system, and it is important to investigate new 
techniques for understanding and quantifying the role of 
entanglement in such systems. Desirable features of these 
techniques include that they be easily computable, even 
for large systems, that they be applicable at finite tem- 
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perature and that they be in principle easy to measure 
and related to known physics. Because most quantum 
systems are described by mixed states, these criteria lead 
naturally to the theory of mixed state entanglement. 

Surprisingly, even the question of whether a mixed 
state of a quantum system is entangled or not is a difhcult 
and much studied question. We refer the reader to the 
many reviews for the literature on the so-called separa- 
bility problem jHHHEil. The difficulty of this problem 
is one of the reasons why computing measures of entan- 
glement can be so challenging and why it is important to 
find more tractable ways of understanding the entangle- 
ment in real physical systems. In this work apply results 
from the theory of the separability problem to study the 
entanglement of quantum many-body systems. This in- 
vestigation leads both to an understanding of the kind of 
Hamiltonians possessing strongly entangled thermal and 
low energy states and also to interesting connections with 
properties of spins systems studied in more conventional 
condensed matter approaches. 

The specific concept that we use from the theory of 
separability is that of an entanglement witness: an ob- 
servable whose expectation value is positive if the state of 
interest is not entangled but for which a negative expec- 
tation value indicates that the state is entangled. An ex- 
ample of such an observable is the Bell observable which 
describes the outcomes of a test for the violation of Bell 
inequalities. In this paper we explore the idea of inter- 
preting Hamiltonians with entangled ground states as en- 
tanglement witnesses. This point of view has attracted 
interest recently. During the preparation of this paper, 
related investigations appeared by Brukner and Vedral 
|40| and by Toth j41,] in which this type of entangle- 
ment witness is studied. As emphasized by Brukner and 
Vedral l40|, because energy is a macroscopic thermody- 
namic property it is reasonable to expect that it could be 
measured in experiment. 

In this paper, we develop the idea of using energy as an 
entanglement witness for quantum many-body systems. 
We introduce two related concepts inspired by the theory 
of entanglement witnesses, and discuss their relevance to 
both ground-state and finite-temperature properties of 
quantum many-body systems. The first is the entangle- 
ment gap: the difference in energy between the ground- 
state energy and the minimum energy that any separable 
(unentangled) state may attain. If the entanglement gap 
of a system is non-zero, then the entanglement of certain 
mixed states is detected simply by measuring their en- 
ergy to be below this threshold. Roughly speaking, if the 
entanglement gap is small then a separable state can be a 
good approximation to the ground state, and we expect 
approximation schemes based on separable states to pro- 
duce reliable results. We investigate how large this gap 
can be for two-spin systems and how this gap depends 
on the coordination number for lattices of coupled spins. 

One advantage of using ideas from studies of mixed- 
state entanglement is that it is possible to investigate 
systems at finite temperature. The second concept we 



introduce is the use of temperature as an indicator of 
entanglement in the thermal state. By comparing the 
thermal energy with the entanglement gap we obtain a 
temperature threshold, the entanglement- gap tempera- 
ture, below which the thermal state is certainly entan- 
gled and we may expect entanglement to influence ther- 
modynamic properties. We show that this temperature 
can become arbitrarily large as the dimension of two in- 
teracting spins increases even if the energy range of the 
system is kept fixed. 

We begin in Sec. with the observation that Hamil- 
tonians with entangled ground states may be viewed 
as entanglement witnesses. We introduce the notion of 
entanglement gap and provide necessary and sufficient 
conditions for this gap to be non-zero. We construct 
a one-to-one mapping between Hamiltonians with non- 
zero entanglement gap and entanglement witnesses. We 
prove a theorem that identifies a set of Hamiltonians with 
the largest possible gap; for spin-1/2 particles, one such 
Hamiltonian is the Heisenberg antiferromagnet. We then 
formally define the entanglement-gap temperature and 
investigate conditions that lead to a high value of this 
temperature. Somewhat counterintuitively, the Hamilto- 
nians with the largest entanglement gap do not have the 
largest entanglement-gap temperature in general. 

In Sec. mil we study the entanglement gap in many- 
body systems, in particular spin models on lattices. We 
prove a general result that the entanglement gap must 
go to zero with increasing coordination number on a bi- 
partite lattice with a fixed local interaction. This result 
is suggestive of a relationship to the success of mean- 
field theory on lattices with high coordination number. 
In Sec. IIVI we conclude by investigating the dependence 
of the entanglement gap on frustration for the Heisen- 
berg antiferromagnet. We show that for such systems 
it is possible to determine that the system is entangled 
even when the reduced state of nearest neighbor spins 
is not entangled. We also investigate how the entangle- 
ment gap behaves near the quantum phase transition in 
the XY model and discuss its relationship to previous 
studies of entanglement at this transition. 

II. HAMILTONIANS AS ENTANGLEMENT 
WITNESSES 

In this section we establish a formal connection be- 
tween Hamiltonians with the property that all low-energy 
states up to a certain energy are entangled, and entan- 
glement witnesses. 

A multipartite mixed state of n subsystems is said to be 
separable if it can be expressed as a convex combination 
of pure product states 

p = E^''iV'')(^'i ^ ^ • ■ • ^ m)m , (1) 

i 

where \ipf) are pure states in the Hilbert space Tij of 
subsystem j, and pi > 0, J^iPi = 1- ^ state can be 
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decomposed in this way, then all correlations are purely 
classical; if not, then there exist truly quantum correla- 
tions and we say that the state is entangled. 

An entanglement witness, Zew, on a multipartite sys- 
tem is a Hermitian operator (observable) with the prop- 
erties that its expectation value in any separable state is 
greater than or equal to zero 

tT[ZEWPscp] > , V pscp € 5 , (2) 

where S is the set of all separable states, and that there 
exists an entangled state, Pcnt, such that 

tr[ZEWPcnt] < . (3) 

We say that Zew witnesses the entanglement of pont- 

For a multipartite Hamiltonian, H, we define the min- 
imum separable energy, 

Escp = min tr[Hpscp] ■ (4) 

Due to the convexity of the set of separable states, this 
minimum can always be achieved by a pure separable 
state. Note that there may be many separable states 
achieving this minimum separable energy i?scp- 

If Escp is strictly greater than the ground-state energy, 
Eq, there is a finite energy range over which all states are 
entangled. We refer to the size of this energy range as 
the entanglement gap. 

Definition: For any multipartite Hamiltonian, H, we 
define the entanglement gap, 

Ge — Escp — Eq , (5) 

where Eq is the ground-state energy of H . The entan- 
glement gap is the energy gap between the ground-state 
energy and the minimum energy that a separable state 
can attain. 

If H has an entangled non-degenerate ground state 
\Eq) then any separable state written in terms of the 
eigenstates of the Hamiltonian must contain contribu- 
tions from higher energy states and must therefore have 
higher energy than Eq . If the ground state is degenerate 
the same argument requires that the entanglement gap is 
greater than zero if there is no state in the ground-state 
manifold that is a pure product state. Conversely a non- 
zero entanglement gap requires that Esep > Eq and so 
there can be no pure product state in the ground-state 
manifold because such a state would have energy Eq. So 
whether or not the entanglement gap is zero depends only 
on the ground-state manifold. A Hamiltonian H has a 
non-zero entanglement gap if and only if no ground state 
of H is separable. 

Constructing Hamiltonians with a non-zero entangle- 
ment gap is straightforward. Every entanglement wit- 
ness can be regarded as a Hamiltonian for a multipartite 
quantum system. For such Hamiltonians, ii^sep = and 
Eq is the minimum eigenvalue of Zew ■ The definition of 



entanglement witnesses implies that Eq < and thus the 
entanglement gap is non-zero. 

Every Hamiltonian with a positive entanglement gap 
Ge > defines an entanglement witness, 

Zew = -ff — Escpl , (6) 

where / is the identity on the total Hilbert space. Be- 
cause Escp is the lowest possible energy for a separable 
state we have tT[ZEwPscp] = tr[i//9scp] — Escp > 0. On 
the other hand if po is a state in the ground-state man- 
ifold we have tr[Z£VKPsep] = Eq — Escp < so Zew is 
an entanglement witness. Note that if H' and H dif- 
fer only by an additive constant they lead to the same 
entanglement witness. We regard such Hamiltonians as 
equivalent. 

In summary, there is a one-to-one map between entan- 
glement witnesses and the equivalence classes of Hamil- 
tonians with non-zero entanglement gap. 

The entanglement gap quantifies the range of energies 
over which all states are necessarily entangled. Note, 
however, that higher energy states may still be entangled. 
So, for example, the thermal state for H must be entan- 
gled for all temperatures such that the thermal energy is 
below Escp but at higher temperatures the thermal state 
may or may not be entangled. 

In appendix^we describe an efficient numerical proce- 
dure, a sequence of semidefinite programs, for evaluating 
the entanglement gap and discuss the concept of bound 
entanglement in this context. 

A. Hamiltonians that Maximize the Entanglement 

Gap 

Having defined the entanglement gap it is natural to 
identify Hamiltonians that have the largest possible en- 
tanglement gap for a given multipartite quantum system. 
We proceed by proving two lemmas, one that the entan- 
glement gap is invariant under local unitary transforma- 
tions of the Hamiltonian, and the other regarding the 
optimal arrangement of the energy levels. We use these 
two lemmas to prove the main theorem of this section, 
which is that a set of Hamiltonians with maximum pos- 
sible entanglement gap are those with a non-degenerate 
maximally entangled ground state and all other eigen- 
states at equal energy. 

Lemma 1. Given a multipartite Hamiltonian, H, and a 
local unitary ?7iocai = C^i <8) C/2 ® • • • ® Un acting on each 
subsystem, the Hamiltonian H' = f^iocai^^^ocai ^"^^ 
same entanglement gap as H . 

Proof: From the cyclic property of the trace, we have 

tr[-H^'Pscp] = ti-[Hpscp], where p^cp = ^LaiPsopC^iocai is 
also separable. That is, for each psep with a certain en- 
ergy under H there is a separable state Pggp with the 
same energy under H' . Therefore Escp = E'^cp- Also, 
because H and H' are related by conjugation by a uni- 
tary they have the same spectrum, and in particular the 
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same ground-state energy. Hence H and H' have equal 
entanglement gap. □ 

We now determine which Hamiltonians have the 
largest entanglement gap. For a comparison of gaps to be 
a sensible, we need to scale by the overall energy range of 
the system. We define the scaled entanglement gap, gs 
as 

gE — Gfi/E^tot , (7) 

where i?tot — -^max ~ Eq is the total energy range, and 
Emax is the highest energy eigenvalue. 

Lemma 2. For any Hamiltonian with scaled entangle- 
ment gap gE, the Hamiltonian H' — I — \Eq){Eq\, where 
\Eq) is a ground state for H , has a scaled entanglement 
gap g'^ greater than or equal to gs- 

Proof: We scale the Hamiltonian so that its lowest eigen- 
value is zero and its highest eigenvalue is one, and thus 
the energy eigenvalues lie in the range < Ei < 1, 
i — . . . dx — where c?t is the dimension of the total 
Hilbert space. The entanglement gap of the scaled 
Hamiltonian H is equal to the scaled entanglement gap 
gE of the original Hamiltonian H. Note that the Hamil- 
tonian H' is already scaled in this manner, i.e., g'^ ~ G'^. 

To prove the lemma, it is sufficient to show the stronger 
result tr[5p] < tr[i7'p], V p, i.e., all states have higher 
energy under H' than under H . To this end 

t-,[pH] = E^mp\E^) < iE^\p\E^) = tr[i?'p] , 

1=0 2=1 

(8) 

as required. The inequality follows from the assumed 
range of energies < Ei < 1 \/ i (where Eq = 0) and the 
fact that < {4>\p\tp) < 1 for any density operator, p, 
and any state, 1-0) (because < p < I)- Therefore E!.^p 
is necessarily greater than Escp (even if the minimum- 
energy separable states are different), and because Eq = 
Eq — and both Hamiltonians are scaled appropriately, 
we have ><?£;, as required. □ 

Using the geometric measure of entanglement for mul- 
tipartite systems defined in '49l| , we consider multipartite 
pure states that are maximally entangled in the sense 
that they have minimal overlap with any pure product 
state, i.e., that they maximize the entanglement measure 

M(|*» = l- max(*|pscp|*). (9) 

Psopes 

Let A/™'^'^ = A/(|^'nic)) be the maximum value of this 
measure, achievable by a maximally entangled state 

|*mc). 

Theorem 1. The largest possible scaled entanglement 
gap of a multipartite system is g'^'^^ = M"^^^ , and can 
he achieved by any Hamiltonian of the form H' ^ I — 
|^mc)(^mc|; where l^'mo) is a maximally entangled state 
by the measure of Eq. 



Proof: The proof follows from the definition of the en- 
tanglement gap, 

gE = l- max (£'o|pscp|-Bo) , (10) 

PaepeS 

and from Lemma |2J □ 

Although we do not present the result here, it is also 
possible to show "so] that all Hamiltonians that have this 
maximum entanglement gap are of this form. 

For multipartite systems it is not known which states 
are maximally entangled according to the measure M . 
However, in |40| examples of highly entangled states are 
given, which place lower bounds on the maximum size 
of the scaled entanglement gap. For example, if each 
of the n subsystems are d-dimensional, there exists a 
symmetrised state \S[n,d)) such that M{\S{n,d))) ap- 
proaches 1 as d"^" in the n ^ oo limit. If each of 
the n subsystems are n-dimensional, there is an antisym- 
metrised state, \A{n)) such that M{\A{n))) = 1 - 1/n!. 
It is clear that entanglement gap can be a very large 
fraction of the total energy range for large numbers of 
coupled systems. 

Bipartite entanglement is much better understood 
than multipartite entanglement, and the following Corol- 
lary gives an explicit form for the maximally entangled 
ground state and the corresponding maximum possible 
scaled entanglement gap for bipartite systems. 

Corollary: The largest scaled entanglement gap for a 
bipartite system Ha ® Hb is gE = ^ — ^/d, where 
d = Toin^dA, ds) is the smaller dimension of the two sub- 
systems, and is achieved by any Hamiltonian of the form 

H' ^ I - \(t)d){(t)d\, where \(j)d) = 7^ Z^iLi Na) and 
are orthonormal bases for Ha/ b- 

Proof: It follows from the convexity of the set of separa- 
ble density matrices that the maximum overlap between 
a pure ground state and a separable state is achieved 
by a pure product state psep = |^)(^| \B){B\, where 
\A) £ Ha, \B) g Hb- In fact, the maximum is achieved 
by setting \A) = \B) = |ls), 

max (Eo\pscp\Eo) ^ Xl , (11) 

PaopeS 

where the Schmidt decomposition 3] for \Eo) is \Eo) = 
X]j=i Xi\iA)\iB), and Ai is the largest Schmidt coefficient. 

Thus, the largest scaled entanglement gap results from 
finding \Eo) with the smallest possible Ai. Normaliza- 
tion iJ2i^i = 1) requires that A^ > 1/d and = l/d 
is achieved by any maximally entangled bipartite state 
|£;q) = \(j)d)- Thus, the Hamiltonian H — I — \4>d){4>d\ 
achieves the maximum possible scaled entanglement gap, 
gE = l- □ 

For dA — dB — 2, the Hamiltonian H — I — \4>2){4'2\, 
where \(j)2) = '^i^i\iA)\iB) / is any maximally en- 
tangled state, has the largest entanglement gap. If 
the Hilbert space corresponds physically to two spin- 1/2 
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systems, then a particularly enlightening example of a 
Hamiltonian of this form is a shifted and scaled version 
of the antiferromagnetic Hamiltonian, H = (Ja-i^b, where 
(Ji, i — A^B \s the vector of Pauli matrices on Tii. It is 
straightforward to show that 



H = I- \ = {aA-cJB + 3/)/4 , 



where \'ip 
state. 



(|0)a|1)s - |1)a|0)b)/\/2 is the singlet 



B. Entanglement-Gap Temperature 

In the following, we investigate the temperature at 
which the thermal state reaches the minimum separa- 
ble energy. We find the temperature below which the 
thermal state is guaranteed to be entangled; this tem- 
perature also provides a non-trivial lower bound on the 
temperature above which the thermal state is guaranteed 
to be separable. 

A quantum system with Hamiltonian, H , in thermal 
equilibrium at temperature, T, is described by the ther- 
mal state 



Pt =exp(-/3i7)/Z, 



(12) 



where f3 — l/kgT is the inverse temperature, kB is Boltz- 
mann's constant, and Z — tr[exp(— is the partition 
function. The energy of the thermal state, the thermal 
energy, is given by 



U{T) = tv[HpT] = - 



1 dZ 
Z'dp 



(13) 



Definition: Given a system with an entanglement gap 
greater than zero, G^; > 0, we define the entanglement- 
gap temperature, Te, to be the temperature at which the 
thermal energy equals the minimum separable energy, 

U{Te) - S,ep. 

The thermal energy is a monotonically decreasing func- 
tion of f3 (i.e., it decreases as the temperature decreases). 
By definition, all states with energy less than Egcp are 
guaranteed to be entangled, and thus the system is cer- 
tainly entangled below the entanglement-gap tempera- 
ture. That is, if we cool our system down below the 
entanglement-gap temperature we know it must be in an 
entangled state. The thermal energy of the system, which 
depends only on the temperature, becomes a witness to 
the entanglement of the thermal state. 

In order to compare Hamiltonians with different to- 
tal energy ranges, -Btot, it is sensible to define a scaled 
temperature, t as 



t ^ keT/Eu 



(14) 



The corresponding scaled entanglement- gap temperature 
is tE = ksTE/Etot- 



For the class of Hamiltonians identified in Theorem 
n with maximal entanglement gap, i.e., H = I — 
|^mc)(^'mc|, where l^'mc) is a maximally entangled state 
by the measure © , it is straightforward to calculate the 
entanglement-gap temperature. The thermal energy is 
given by 



Uit) 



(rfT-l)exp(-/3) 
l + {dT- 1) exp(-/3) ■ 



Setting U{tE) = Esep = M'"^'' gives 



tE 



{dr - 1)(1 - M'"^'^) 



(15) 



(16) 



As an example, we consider the entanglement-gap tem- 
perature of a bipartite system (each subsystem of dimen- 
sion d), with Hamiltonian H = /— and scaled en- 
tanglement gap = 1 — 1/d. The scaled entanglement- 
gap temperature for this system is 



i£ = [log,(d+l)] ' 



(17) 



Note that the entanglement-gap temperature decreases 
with increasing dimension despite the fact that the en- 
tanglement gap increases. This behavior is due to the 
fact that the number of eigenstates with energy one in- 
creases quadratically with d, while the ground state re- 
mains non-degenerate. 



C. Hamiltonians of Bipartite Systems Possessing 
Large Entanglement-Gap Temperature 

It is natural to ask which Hamiltonians have the high- 
est (scaled) entanglement- gap temperature. Somewhat 
counterintuitively it is not the Hamiltonians with the 
largest entanglement gap. In fact, there are Hamiltonians 
with arbitrarily high entanglement-gap temperature. To 
provide an example, we restrict our attention to the case 
where the two subsystems of the bipartite system have 
the same dimension, dA = ds = d. The projectors onto 
the symmetric and antisymmetric subspaces of Ti.A ^Ti-B 
areHs = {I-\-V(^a,b))/2, andH^ = (/- V(a,b))/2, respec- 
tively, where V(a,b) is the permutation operator on the 
two subsystems, defined by V(a,b)\^)a\4>)b = \4>)A\i')B 
for all 10). The antisymmetric subspace contains only 
entangled states. Thus if we define a Hamiltonian as the 
projector onto the symmetric subspace 



H = ns, 



(18) 



then all symmetric states have energy one, all antisym- 
metric states have energy zero, and there is a finite en- 
tanglement gap. We find the gap by directly calculating 
the energy of a pure separable state, |yl)|i?). 



{A\{B\H\A)\B) = il + \{A\B)n/2. 



(19) 



From this expression it is clear that the minimum en- 
ergy of 1/2 is achieved by any pure separable state such 
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that {A\B) = 0. The entanglement gap is Ge = 1/2, 
independent of d. 

For the symmetric-projector Hamiltonian the thermal 
energy is given by 



m = 



d{d + 1) exp(-/3) 



d{d~ l) + d{d+l)exp{~(3) 
Using i^sop = 1/2 we find 



d/2, for d > 1 . 



d+1 



-1 



(20) 



(21) 



Remarkably, for this Hamiltonian the scaled 
entanglement-gap temperature increases without 
bound as the dimension of the subsystems increases. 

Thus, for Hamiltonians that only have eigenvalues 
zero or one, there is a trade-off between ground-state 
degeneracy and the entanglement gap in determin- 
ing the entanglement-gap temperature. Even though 
the Hamiltonian with a non-degenerate maximally- 
entangled ground state has a larger entanglement 
gap, the symmetric-projector Hamiltonian has a higher 
entanglement-gap temperature due to its large ground- 
state degeneracy. 

In Appendix^ we investigate other Hamiltonians with 
ground-state manifolds containing only entangled states, 
and present evidence that no other bipartite Hamil- 
tonian with a two-level energy spectrum possesses an 
entanglement-gap temperature greater than the Hamil- 
tonian (|18|l . In Appendix |0 we investigate the entan- 
glement temperature of two qubit systems, and provide 
evidence that the Heisenberg antiferromagnetic Hamilto- 
nian has the highest entanglement-gap temperature. 

We note that Toth gives an example of a multi- 
party Hamiltonian, the Heisenberg interaction between 
all pairs of n spin- 1/2 particles, whose entanglement-gap 
temperature increases linearly with n, i.e. it is arbitrar- 
ily high for arbitrarily large systems. However, unlike 
our example, the total energy range also increases lin- 
early with n. The scaled entanglement-gap temperature 
of their Hamiltonian therefore approaches a constant as 
n —^ oo. By contrast, the entanglement-gap temperature 
of our example is arbitrarily high despite the fact that 
the total energy range is bounded. 



III. THE ENTANGLEMENT GAP OF 
QUANTUM MANY-BODY SYSTEMS 

In this section, we investigate the entanglement gap 
for quantum systems arranged on some graph or lattice 
that interact with some local neighborhood. Because we 
are only considering finite-dimensional systems, the sub- 
systems can always be thought of as spins of some total 
angular momentum, so we use the terms "subsystem" 
and "spin" interchangeably. For a particular type of cou- 
pling - bipartite lattices - we provide an explicit method 



for calculating the entanglement gap, which applies to 
various spin systems often considered in the condensed 
matter literature. We also prove that, as the coordina- 
tion number grows, the entanglement gap per interaction 
must decrease to zero. This result makes use of the fact 
that as the number of equivalent spins connected to a 
given spin in the lattice increases there does not exist a 
global state of the system for which each interacting pair 
is strongly entangled. 



A. The Entanglement Gap of 2-Local Hamiltonians 

We now consider multipartite systems with only two- 
body interactions. The Hamiltonian for such a system 
can be defined by a set of coupling Hamiltonians Hij that 
act as the identity on all the spins other than i and j, 
and a graph or lattice where the vertices represent spins 
and edges represent an interaction between the spins on 
the two sites. We refer to each two-body interaction, or 
edge on the graph, as a bond. The Hamiltonian for the 
entire system is 



H = 



E 



(22) 



where J2<ij> indicates a sum over vertices connected 
by an edge, i.e., a sum over bonds. We refer to such a 
Hamiltonian as 2-local. It follows that that the energy 
of a 2-local Hamiltonian depends only on the reduced 
density matrices of each interacting pair; see e.g. |10|. 
Thus, the energy for any global state p is 



E = ti[Hp] = tr[i?»jPij] 



(23) 



<ij> 



where is the reduced state of the interacting pair of 
spins <i,j>. In the following wc consider only systems 
where each of the coupling Hamiltonians Hij is equal. 

We note that the reduced states pij are not completely 
arbitrary if they are to be consistent with a global state 
p for the whole system. In particular, the ground state of 
the graph or lattice cannot simply be constructed from 
the reduced states that minimize the energy of each bond 
{tr[Hij Pij]) , because these reduced states may not be con- 
sistent with a global state. As we demonstrate below, 
this situation can occur when there is a non-zero entan- 
glement gap for the coupling Hamiltonian Hij. Moti- 
vated by the results of Sec. ^1 well as its importance 
in condensed matter physics, we use the Heisenberg an- 
tiferromagnet as our standard example. 

A bipartite graph or lattice is one for which the ver- 
tices can be divided into two sets, A and B, such that 
the edges only connect vertices from A with vertices from 
B. Examples of bipartite lattices include the square lat- 
tice, see Fig. Q] and the hexagonal lattice on the plane. 
An even number of vertices arranged in a ring is also 
bipartite. 
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For bipartite graphs or lattices, we now demonstrate 
how to construct a separable state with the lowest pos- 
sible energy. First, consider a minimum-energy sep- 
arable state for a pair of spins under the in- 
teraction Hamiltonian, and construct a global state, 
^iAeA\A)ij^ ^i^eB \B)i^, such that all the spins on the 
subset A are in the state \A) and likewise for B. By 
Eq. H23|) the energy per bond of this state is the same 
as the energy E'sep of the state and it provides 
an upper bound on the minimum separable energy per 
bond of the full Hamiltonian H22|l. To see that there 
is no global separable state with lower energy, suppose 
that such a state exists. Then all of the nearest neighbor 
reduced density matrices pij must be separable and by 
Eq. I|23|) at least one of them must have a lower energy 
under the interaction Hamiltonian than a con- 
tradiction. Therefore, the state ®iBeB \B)iB 
is indeed a minimum-energy separable state of the entire 
system. 




FIG. 1: Examples of n-partite lattices, (a) square lattice 
(n = 2), (b) triangular lattice (n — 3), (c) kagome lattice 
(n = 3), (d) checkerboard lattice (n = 4). The n differ- 
ent markers indicate the n subsets that the vertices of the 
n-partite lattice may be divided into so that there are only 
interactions between distinct subsets. 



Given this result, we can determine the lowest possible 
energy of a separable state for the Hamiltonian (|22(l on 
any bipartite graph or lattice simply by solving the prob- 
lem for a single pair of spins. Finding the entanglement 
gap for such systems reduces to finding the entanglement 
gap of the two-body interaction and the ground-state en- 
ergy Eq of the overall system. This important fact was 
noted by Toth i41] who provided a slightly different ar- 
gument. 

Consider a bipartite graph with Hamiltonian H. Then 
our argument proves that the operator H — NEgcp is an 
entanglement witness, where N is the number of bonds. 
Note that we can express this entanglement witness as a 



sum over bonds, 

H - NE,,p = J2 - ^-p) ' (24) 

<i,J> 

where each term in the sum is a bipartite entangle- 
ment witness. As a result the expectation value depends 
only on the bipartite reduced density matrices of nearest 
neighbors and can only be negative if these reduced den- 
sity matrices are entangled. This result can be extended 
to apply to lattices (with N oo) as well. So, while 
the ground-state energy is certainly a global quantity, 
this construction is only sufficient to detect the existence 
of bipartite entanglement between interacting pairs in a 
bond. 

In a similar way we can calculate the entanglement gap 
for n— partite graphs or lattices that are formed by groups 
of n spins each having an "all-to-all" interaction graph. 
We first construct a minimum-energy separable state of 
a single group of n spins with an all-to-all interaction 
graph. This state extends to minimum-energy separa- 
ble state of the entire lattice as depicted in Fig. We 
then compare the minimum-separable energy per bound 
to the ground state energy per bond to find the entan- 
glement gap per bond. For example, on the tripartite 
triangular and kagome lattices (see Fig. it is possible 
to calculate the entanglement gap from the ground-state 
energy and the minimum separable energy of three spins 
having an all-to-all interaction graph (a single triangle). 
We describe the nature of the entanglement that can be 
witnessed in these systems in Sec. IIVBI 



B. Entanglement Gap and coordination Number 

The coordination number of a lattice is the number of 
edges incident on each lattice site, i.e., the number of 
other systems that each spin interacts with via the cou- 
pling Hamiltonian Hij. As we are now considering lat- 
tices, our assumption that all interactions Hij are equal 
implies translational symmetry. We now investigate how 
the entanglement gap varies with the coordination num- 
ber of the lattice. The basic idea stems from the fact that, 
as a result of the translational symmetry, the ground 
state of our 2-local Hamiltonians have equal reduced den- 
sity matrices for interacting pairs ^ . As the coordination 
number increases, this equality requires that every spin 
shares the same reduced density matrix with an increas- 
ing number of other spins. The results of |^ [s^ 
then preclude the reduced density matrices from being 
strongly entangled. Building on these results, we prove a 
theorem stating that, as the coordination number of the 
lattice grows, the entanglement gap decreases to zero. 



If there is spontaneous symmetry breaking a mixture of the 
symmetry-broken ground states will be translationally invariant. 
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We then investigate this behaviour in the specific case of 
the Heisenberg antiferromagnet. 

In order to prove resuhs about the maximum possible 
entanglement gap in Sec. Ill Al it was natural to use the 
scaled entanglement gap. However in what follows it is 
more convenient to use the entanglement gap per bond, 
Ge/N, where N is the total number of bonds. (Note that 
this entanglement gap per bond is well defined even for 
lattices with TV — > oo.) These two methods of scaling are 
roughly equivalent because the total energy range of the 
system tends to scale linearly with the number of sites. 

We begin by considering a restricted set of graphs 
which we will use to prove results that bound the en- 
tanglement gap on any bipartite lattice. We define a star 
graph as a bipartite graph where there is only a single ver- 
tex, the centre, in one subset, A = {^o}i a-nd k vertices, 
the points, in the other subset, B — {Bi, i — 0, . . . , k—1}, 
and where edges connect each point and the centre. 

A strictly positive entanglement witness, Zewi is a- 
Hermitian operator whose average is strictly positive on 
separable states, tT[ZEwPsep] > 0, V p^cp G S, but 
which has at least one negative eigenvalue. 

Before stating and proving our main theorem we 
present the following lemma: 

Lemma 3. Let Zew &e a stictly positive entanglement 
witness acting on TLao ®T^Bo ■ Then there exists a positive 
integer k such that 

E ViBoMZEW ®tl lB,)yiBo,B.) > , (25) 

4=0 

where V(Bo.Bi) the self-adjoint unitary operator that 
swaps the Hilbert spaces 'Hbo cind Tis- . 

This Lemma is a straightforward modification of The- 
orem 2 of j45i and the proof follows similarly ^. 

Using the general mapping between entanglement wit- 
nesses and Hamiltonians with non-zero entanglement gap 
discussed in Sec. this result on strictly positive en- 
tanglement witnesses bounds the entanglement gap for 
Hamiltonians on star graphs. 

Theorem 2. For any coupling Hamiltonian HaqBo <^'^d 
any e > there exists a positive integer, k, such that 
the entanglement gap per interaction for the Hamiltonian 
\2^) on a star graph with k points is less than e. 

Proof: The non-trivial case occurs when HaqBq has non- 
zero entanglement gap. Note that the total Hamiltonian 
on the star graph may be written as 

fc-i 

TJ^t- = E V(Bo,B.) {HaoB„ ®tl Ib,)V;^^^s,) ■ (26) 

i=0 



^ The key difference between this proof and the one found in |45] is 
that here, mixing with the operators Vj^g forces the result- 
ing entanglement witness to be block-diagonal in the irreducible 
representations of the symmetric group, rather than proje ctin g 
into the fully symmetric representation of this group as in l45l . 



We define a strictly positive entanglement witness on 
Ti-Ao T^Bo as 

Zew = HaqBq — Escpl + , (27) 

where E'sop is the energy of the minimum-energy separa- 
ble state of HaqBo and, by adding e > 0, Zew is guaran- 
teed to be a strictly positive entanglement witness. From 
Lemma 121 there exists a k such that 

fe-i 

E ^(So,B.) {{HaoBo - i?sop + e) ®tl Ib.) > > 

i=0 

and so H^^^'^ > k{Esep — e). Because the energy of the 
minimum-energy separable state of H^^'^'^ is kEgcp, this 
implies that the entanglement gap of the total Hamilto- 
nian satisfies G^'^'' < ke. Thus, given any e > there 
exists a k such that G^e^'^ jk < e, as claimed. □ 

As an illustration of this theorem we consider the spin- 
1/2 Heisenberg antiferromagnetic Hamiltonian on a star 
graph. Recall that the coupling Hamiltonian is Hij — 
(fi-dj; the entanglement gap of this coupling Hamiltonian 
was investigated in Sec. IH Al The ground state is the 
singlet, with energy —3, and the three triplet states all 
have energy -1-1. The minimum-energy separable states 
are of the form \ A) \B) such that {A\B) — 0, with energy 
-1. 

Using the permutation symmetry amongst the points 
of the Hamiltonian on the star graph it is possible to cal- 
culate its ground-state energy exactly 53], Eq — —(k+2). 
The coordination number of the center of the graph is the 
number of points, k. The energy of any minimum-energy 
separable state is Escp = —k. In Table we present the 
entanglement gap per bond and the scaled entanglement 
gap for comparison with other lattices below. 



Coord. 


Eo per 


Esep per 


Ent. Gap 


Scaled 


No. k 


bond 


bond 


per bond 


Ent. Gap 


1 


-3 


-1 


2 


0.5 


2 


-2 


-1 


1 


0.333 


3 


-1.667 


-1 


0.667 


0.25 


4 


-1.5 


-1 


0.5 


0.2 


5 


-1.4 


-1 


0.4 


0.167 


6 


-1.333 


-1 


0.333 


0.143 



TABLE L Properties of star graphs with the Heisenberg anti- 
ferromagnetic Hamiltonian as a function of coordination num- 
ber k. The ground-state energy, minimum separable energy 
and entanglement gap are all per bond, i.e., energies divided 
by k. The scaled entanglement gap is the entanglement gap 
divided by the total energy range of the system. 

Although the Heisenberg antiferromagnet has the 
largest entanglement gap for two qubits, we have not 
proved that it has the largest entanglement gap per bond 
on a star graph. However we have calculated the entan- 
glement gap per bond for numerous common spin mod- 
els such as the XX Z model and XY model, all of which 
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have a smaller entanglement gap per bond. Therefore we 
conjecture that the Heisenberg antiferromagnet has the 
largest entanglement gap per bond on a star graph. If 
this were true it would provide an upper bound of 0(1 /k) 
on the approach to zero of the entanglement gap per bond 
implied by Theorem [3 

In order to determine the entanglement gap on a bi- 
partite lattice, we require the ground-state energy of the 
lattice as well as the lowest energy achievable by a sep- 
arable state of a single pair of spins, as noted above. 
The ground-state energy of a star graph can be used to 
bound the ground-state energy of a bipartite lattice with 
the same coordination number as follows. 

Lemma 4. The ground-state energy per bond of any cou- 
pling Hamiltonian on a bipartite lattice with coordination 
number k is greater than or equal to the ground-state en- 
ergy per bond of the same Hamiltonian on a star graph 
with k points. 

Proof: The essential idea is to divide the expression 
for ground-state energy on the lattice into a sum over 
star graphs with k points where k is the coordination 
number of the lattice. Let po denote the translationally- 
invariant ground state of the entire lattice. Consider the 
star graph consisting of a particular lattice site (the cen- 
ter) and those sites connected to it by a coupling (the 
points). The reduced state on the star graph is obtained 
by tracing out all sites not in the star 

Pstar = trj^lstar} [Po] ■ (28) 

This state is independent of the lattice site chosen as the 
center (due to the translational invariance of po), and 
the energy per bond of this reduced state is the same 
as the ground-state energy per bond of the lattice. The 
ground-state energy is then 

Eo = ti[Hpo] = J2 i^Hste^rPst^A/k , (29) 

i 

Furthermore, the energy of pstar can only be greater than 
the energy of a ground state |i?o)star of the star Hamil- 
tonian, 

tr[77starPstar] > tr[i?star |£^o)star (-Eq |] . (30) 

It follows that the ground-state energy per bond of the 
bipartite lattice is greater than or equal to the ground- 
state energy per bond of the star graph. □ 

We note that a similar argument is used in |54| to 
bound the ground-state energy of the Heisenberg antifer- 
romagnet. 

Using this bound for the ground-state energy, it is 
straightforward to bound the entanglement gap on bi- 
partite lattices, which is the main result of this section. 

Theorem 3. Given any e > there exists a positive 
integer, k, such that the entanglement gap per bond for 
an arbitrary coupling Hamiltonian on any bipartite lattice 
with coordination number k is less than e. 



Proof: Because the bipartite lattice and star graph are 
both bipartite, they have the same minimum separable 
energy per bond. The result now follows from Theorem 
Inland Lemma El □ 



Lattice 


Coord. 


Eo per 


iJscp per 


hAa. Cjap 


scaled 




iNo. 


bond 


bond 


per bond 


Ent. Gap 


single bond 


1 


-3 


-1 


2 


0.5 


ID chain 


2 


-1.772 


-1 


0.772 


0.279 


hexagonal 


3 


-1.452 


-1 


0.452 


0.184 


square 


4 


-1.338 


-1 


0.338 


0.145 


cubic 


6 


-1.194* 


-1 


0.194 


0.088 


single triangle 


2 


-1 


-0.5 


0.5 


0.25 


kagome 


4 


-0.874 


-0.5 


0.374 


0.200 


triangular 


6 


-0.726 


-0.5 


0.226 


0.131 


single 
tetrahedron 


3 


-1 


-0.333 


0.667 


0.333 


checkerboard 


6 


-0.67^ 


-0.333 


0.34 


0.20 



TABLE II: Entanglement gap for the Heisenberg antiferro- 
magnet for various bipartite and frustrated lattices with dif- 
ferent coordination numbers. Ground-state ener gies taken 
from JSa], except * from linear spin- wave theory and ^ 
from ISTII. 



To illustrate this theorem, we calculate the entangle- 
ment gap per bond of the spin-1/2 Heisenberg antiferro- 
magnet on simple bipartite lattices with varying coordi- 
nation number. In Table Ull we present the ground-state 
energy, taken from the literature, and thus the entangle- 
ment gap per bond and scaled entanglement gap for a 
Heisenberg antiferromagnet on a ID chain, honeycomb, 
square and cubic lattice (all bipartite), as well as some 
non-bipartite lattices to be discussed in Sec. IIVBI It 
can be seen that the entanglement gap per bond does 
decrease with increasing coordination number for the bi- 
partite case, and is always less than that of the corre- 
sponding star graph in Tabled as proved by Lemma 0] 
The entanglement gap per bond appears to decrease with 
coordination number on tripartite lattices as well, thus 
providing evidence that this behavior is not confined to 
bipartite lattices. 

IV. DISCUSSION 

In this section, we discuss some of the implications of 
our results and explore the connections with other results 
from the condensed matter literature. We also discuss 
frustrated lattices and quantum phase transitions and 
their effect on the entanglement gap. 

The energy gap between the ground-state energy and 
the lowest energy achieved by a separable state has been 
discussed in the quantum magnetism literature using 
a slightly different terminology 55]. There, separable 
states are associated with "classical configurations," ar- 
rangements of classical spin vectors minimizing the en- 
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ergy of the appropriate classical Heisenberg spin model. 
The reduction in ground-state energy below this point 
is typically ascribed to "quantum fluctuations." As a 
result, Table HTl is essentially drawn directly from the re- 
view by Lhuillier and Misguich ;55J. Our results show 
that, in this context at least, the term "quantum fluctua- 
tions" as discussed in the condensed matter literature can 
be identified precisely with entanglement as discussed in 
the quantum information literature, and the associated 
reduction of ground-state energy in antiferromagnets can 
be directl y re lated to the theory of mixed state entangle- 
ment mililal. 

Note that the entanglement gap is well over a quarter of 
the total energy range for Heisenberg antiferromagnet on 
a line and is a very significant fraction of the total energy 
range for the majority of the lattices considered. This 
large entanglement gap reinforces the argument made by 
Brukner and Vedral that the entanglement witnesses 
resulting from the energy of appropriate spin models can 
have macroscopic expectation values. 



A. Bipartite Lattices and Mean-Field Theory 

Mean-field theory is a term used to describe a variety 
of techniques in condensed matter physics for finding an 
approximation to the ground state of a quantum many- 
body system. Typically such techniques correspond to 
searching for a separable state that approximates the 
ground state. It is a well-known observation that mean- 
field theory is more accurate in higher dimensions and, 
because coordination number typically increases with the 
dimension, for higher coordination number. So for exam- 
ple, dynamical mean field theory for fermion systems is 
known to be exact in infinite dimensions |58l |. 

In the context of our present work, we expect mean- 
field theory to work well when the entanglement gap is 
small, because there exists a separable state that has en- 
ergy close to the ground-state energy and thus a varia- 
tional approach involving separable states might be ex- 
pected to be accurate. Theorem |31 demonstrates in a pre- 
cise way that the entanglement gap decreases to an arbi- 
trarily small value with increasing coordination number 
on bipartite lattices, independent of the particular cou- 
pling Hamiltonian. This result is therefore suggestive of 
a quantitative connection between entanglement and the 
improvement of mean-field theory with dimension. 

The work of Raggio and Werner [sj aimed to develop 
a rigorous mean field theory for Hamiltonian models on 
star graphs with a large number of points. Our results are 
ultimately based on a characterization of bipartite sepa- 
rable states proven there and in Ref. '5^, subsequently 
used in Ref. to prove a result closely related to our 
LemmalHl The proofs in Ref. [sj are technically very dif- 
ficult, because they apply not only to finite dimensional 
spin systems but to any quantum system defined on a sep- 
arable Hilbert space. These results may provide a more 
direct route to our Theorem [21 for star-shaped graphs. 



which could then be used to prove the result for bipartite 
lattices in more generality; however, we have preferred to 
give a simple derivation valid for finite-dimensional spin 
systems. 

B. Frustrated Lattices and Multipartite 
Entanglement 

Lattices that are not bipartite lead to spin systems that 
are often referred to as frustrated in condensed matter 
physics [13. This terminology arises from the fact that 
the minimum-energy separable state for two neighboring 
sites on such a lattice is not equal to the minimum-energy 
separable for the two sites coupled alone. As a result the 
energy per bond on such a lattice is higher than the en- 
ergy of a single pair for the same interaction The 
physics of frustrated quantum and classical spin systems 
have been a subject of intensive research in recent years; 
we refer the reader to |59j | for a review. In the follow- 
ing, we briefiy consider the effect of frustration on the 
entanglement gap. 

Further motivation for studying lattices that are not 
bipartite comes from considering the nature of the en- 
tanglement detected by the Hamiltonian. On bipartite 
lattices, entanglement is only detected when the reduced 
density matrices associated with each bond are entan- 
gled. So, for example, states which are multipartite en- 
tangled but contain no bipartite entanglement, such as 
the three-party GHZ state, will never have lower energy 
than the minimum separable energy on a bipartite graph 
for any interaction Hamiltonian. On non-bipartite lat- 
tices it is sometimes possible for a coupling Hamiltonian 
to witness the entanglement of such states. 

A simple example of a non-bipartite lattice is the reg- 
ular triangular lattice, which is tripartite but not bi- 
partite. We consider two other non-bipartite lattices in 
two dimensions: the kagome lattice, which is made up 
of corner-sharing triangles, and the checkerboard lattice. 
These lattices are depicted in Fig. ^ 

Again we consider the Heisenberg interaction. In order 
to find the lowest separable energy for the triangular and 
kagome lattices, we first find a minimum-energy separa- 
ble state for a single triangle, as described in Sec. IIII Al 
The total Hamiltonian for a single triangle is 

iJ = (7l • (72 + <T2 • (73 + (?3 • <Tl . (31) 

Its spectrum and minimum-energy separable states may 
be found by standard symmetry methods (for exam- 
ple 113). The ground state is four-fold degenerate with 
energy Eg — —3. For the Heisenberg antiferromagnet of 
n spins with an all-to-all coupling, it is straightforward to 



^ Dawson and Nielsen l27f derive a bound on the ground-state en- 
tanglement based on the frustration of the quantum Hamiltonian, 
not the frustration of its classical counterpart. 
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show that a minimum-energy separable state is given by 
any configuration of spins where the total spin vector is 
zero [Hflj l. For the triangle, a minimum-energy separable 
state is 

|T)i ® (|T>2 + \/3U)2)/2® (|T)3 - \/3U>3)/2 , (32) 

which corresponds to a classical configuration of spins 
at an angle of 27r/3 from each other in the plane hav- 
ing a total spin zero (the "Mercedes star" configuration 
in li^)- This state has energy -Esep. = —3/2. The max- 
imum energy manifold is spanned by the states with all 
spins parallel, and has energy -Emax = 3. 

From these results we can calculate the entanglement 
gap per bond and the scaled entanglement gap for the 
Heisenberg interaction on the triangle, shown in Table ITll 
Also shown are the entanglement gaps for the kagome and 
triangular lattices, calculated from Egcp for the triangle, 
and the ground-state energy of the entire lattice, taken 
from 55]. Note that, as for bipartite lattices, the entan- 
glement gap per bond appears to decrease with coordina- 
tion number. We have also considered the checkerboard 
lattice (see Fig. ^ which is made up of corner-sharing 
tetrahedra and has a coordination number of six. We ob- 
tain the ground-state energy of this model from Ref . [s^l , 
where it is estimated from exact diagonalization of small 
samples. 

The reduced density matrices associated with bonds 
of the lattice in the ground state are not entangled for 
these frustrated systems. Note that the symmetries of 
the Heisenberg model guarantee that these bipartite re- 
duced density matrices are so-called Werner states [gi) . 
invariant under any local unitary rotation of the form 
U ® U . These states are entirely characterized by the 
fraction of the population that is in the singlet state, and 
when this fraction is less than a half the state is separa- 
ble ■ Therefore the reduced density matrix associated 
with each bond is separable whenever the ground-state 
energy per bond is above the minimum separable energy 
of the Heisenberg Hamiltonian for a single pair of spins. 
With the ground-state energy per bond from Table ITTl it 
is clear that there is no bipartite entanglement of nearest 
neighbor spins for the Heisenberg model on the trian- 
gular, kagome or checkerboard lattices (because per 
bond is greater than -1). The entanglement gap for these 
systems is associated with the entanglement of the re- 
duced states of the triangles or tetrahedrons that make 
up the lattice. Thus, the Hamiltonian serves as a witness 
for multipartite entanglement in these systems. 

It appears that as the frustration of the classical spin 
model increases, so does the entanglement gap. For a 
coordination number of six the entanglement gap as a 
fraction of the overall energy range of the Hamiltonian 
increases from 0.088 on the bipartite cubic lattice to 0.131 
on the tripartite triangular lattice and finally to around 
0.2 on the checkerboard lattice. It would be interesting 
to understand this behavior in more detail. It is a fea- 
ture of frustrated classical spin models that they have a 
large number of configurations achieving the lowest pos- 



sible energy, which may be a contributing factor to this 
observed larger entanglement gap. 



The Entanglement Gap in a Simple Quantum 
Phase Transition 



The role of entanglement in quantum phase transitions 
[ll| is currently of considerable interest P, |E IE ^ • Per- 
haps the simplest model to exhibit a quantum phase tran- 
sition, used in many of these studies, is the ID infinite- 
lattice transverse field XY model with Hamiltonian 



H 



N-l 

E 



1 + 7 



(33) 



where 7 is the anisotropy in the x ~ y plane, and A is 
an external magnetic field, N is the total number of lat- 
tice sites, and cyclic boundary conditions are imposed so 
that a subscript N is identified with 0. For 7 = 1 the 
transverse field Ising model is recovered. 

It is of interest to see how this phase transition affects 
the entanglement gap. Here we calculate the entangle- 
ment gap of the ID XY model as a function of (7, A) 
in the thermodynamic {N — > 00) limit. Because a \D 
lattice is bipartite (for N even), given knowledge of the 
ground-state energy it is sufficient to calculate the en- 
tanglement gap for the coupling Hamiltonian in order to 
calculate the entanglement gap of the entire system, as 
described in Sec. IIIII In this case the coupling Hamilto- 
nian may be chosen to be 



TjXY 1+7 X^X 



where the factor of 1/2 in front of the magnetic field 
accounts for the fact that each site is involved in two 
interactions. In Appendix^ we calculate the minimum- 
separable energy for this coupling Hamiltonian, Eq. (|D4(I . 

The XY model on a \D chain, Eq. H33|) , is well-known 
to be exactly solvable via the Jordan- Wigner transfor- 
mation; see e.g. ^ij- We obtain the ground-state energy 
from this method. 

In Fig.|21we plot the scaled entanglement gap as a func- 
tion of (7, A) in the thermodynamic limit. The quantum 
phase transition in this model occurs at A = 1 for 7 7^ 0. 
Previous studies have indicated that the ground state be- 
comes highly entangled at this point, and this behavior is 
manifest in a sudden rise in the entanglement gap about 
this point. Intuitively one might expect that the more 
entangled the ground state, the larger the entanglement 
gap. While qualitatively true, this connection cannot be 
exact because the entanglement gap is a property of the 
whole Hamiltonian; it can depend on all energy eigen- 
states and their energies and is not just a property of the 
ground state. 

Nevertheless, given the discussion in Sec. Ill Al it is rea- 
sonable to consider a connection between the ground- 
state entanglement as measured by M(|5')) of Eq. 
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FIG. 2: Entanglement gap as a function of anisotropy 7 and 
transverse field A for XY Hamiltonian on a ID lattice in the 
thermodynamic limit. 

and the entanglement gap. The entanglement of the 
ground state under this measure has been investigated 
in the XY model in recent work by Wei et al. This 
measure depends only on the maximum overlap of the 
entangled state with a separable state. One might 
expect that the minimum-energy separable state is one 
which has the maximum overlap with the ground state. 
However, in order to achieve maximum overlap with the 
ground state it may also be necessary for a separable 
state to have large overlap with high-energy eigenstates. 
Therefore a separable state may achieve lower energy by 
having less overlap with the ground state but consider- 
ably more overlap with low-lying excited states. 

In Ref. the derivative of the global entanglement 
with respect to the external field was found to contain a 
singularity at the critical point consistent with the uni- 
versality class of the model. Although we see a qualita- 
tively similar peak near the critical point, there is no sin- 
gularity in the derivative of the entanglement gap. Again, 
such a singularity may not have been expected because 
the entanglement gap is not simply a property of the 
ground state. 

D. Summary 

We have studied entanglement in quantum many-body 
systems from the point of view of the Hamiltonian as an 
entanglement witness. We introduced two related con- 
cepts useful in studying the role of entanglement in the 
ground and thermal states of multipartite quantum sys- 
tems. The first is the entanglement gap, which is the dif- 
ference in energy between the ground-state energy and 
the minimum energy that any separable state can attain. 
If the energy of the system lies within the entanglement 
gap range, the state of the system is guaranteed to be 
entangled. The second concept is the entanglement-gap 



temperature: the temperature at which the energy of the 
thermal state is equal to the minimum separable energy, 
and below which the thermal state must be entangled. 
The entanglement- gap temperature provides a threshold 
for deducing the thermal state of the system to be entan- 
gled, based on its energy. 

For multipartite, finite-dimensional quantum systems 
we proved that Hamiltonians possessing a non-degenerate 
maximally entangled ground state (according to a global 
measure of entanglement) and all other energy eigen- 
states degenerate maximize the entanglement gap. The 
related question of which Hamiltonians have the high- 
est entanglement gap temperature is more challenging; 
substantial evidence is given that the Heisenberg anti- 
ferromagnetic Hamiltonian has the largest entanglement 
temperature for two qubits. 

On bipartite lattices, i.e., those lattices for which there 
are only interactions between two disjoint subsets of the 
vertices, we proved that the entanglement gap decreases 
to zero as the coordination number increases. This re- 
sult suggests a quantitative reason why approximation 
schemes based on separable states, such as various forms 
of mean-field theory, appear to give more reliable results 
at higher coordination number. 

On frustrated lattices, i.e., those that are not bipartite, 
we noted that the Hamiltonian can act as an entangle- 
ment witness for multipartite entanglement, even when 
there is no bipartite entanglement present. Finally, we 
calculated the entanglement gap near a simple quantum 
phase transition, and showed that although it does not 
follow any universal scaling law, it does increase near 
the quantum phase transition, as may have been ex- 
pected from previous studies in which the ground state 
was found to become highly entangled at that point. 
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APPENDIX A: SEMIDEFINITE PROGRAMS 
FOR THE ENTANGLEMENT GAP 

We will now describe efficient numerical procedures for 
evaluating the entanglement gap of a given Hamiltonian 
using semidefinite programs. 

Semidefinite programs are a type of convex optimiza- 
tion problem '4^, H§\ , which are appealing because they 
have efficient numerical implementations. With the view 
of Hamiltonians as ent ang lement witnesses, and following 
methods described in [iSlisllI^ it is possible to express 
the problem of finding the minimum separable energy 
as a sequence of semidefinite programs, whose solutions 
converge to Egep- The simplest program, which applies 
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for bipartite systems with Hilbert space Ha ^H-b, is 

max e , 
subject to H -fl = P + Q^^ , 
P > 0, 

g>0, (Al) 

where Ta denotes the partial transpose over system A. 
Let di = dim(Hi), i = A,B. When = 2 and = 2 
or 3 the maximum e obtained from this program cor- 
responds to the minimum separable energy, i?sop- The 
optimum value e* of the semidefinite program gives an 
entanglement witness Zew = H—e*I, and a lower bound 
on the entanglement gap equal to the largest magnitude 
negative eigenvalue of Zew- The entanglement witness 
produced by (|A1|I is referred to as decomposable because 
it can be written Zew = H-e*I ^ P + Q'^^ for P > 0, 
Q > 0, and can only detect entangled states with non- 
positive partial transpose. 

If the subsystems are of higher dimension, it is pos- 
sible for an entangled state to have a positive partial 
transpose. Such states are hound entangled |47l |. and the 
semidefinite program ljAl|l only finds the gap between 
the ground-state energy and the minimum-energy pos- 
itive partial transpose state. This solution provides a 
lower bound on i?sep, and it is possible to devise a nested 
sequence of programs that provide increasingly tighter 
bounds 0- 

As all entanglement witnesses may be viewed as Hamil- 
tonians with entangled low energy states, one way of pro- 
ducing bound entangled states suggests itself: as thermal 
states. An example of a Hamiltonian for which there are 
bound entangled states which achieve lower energy than 
any separable state may be derived from the Choi form, 
as described in [45ll4^ . This Hamiltonian, which acts on 
the minimal-dimension system on which bound entangled 
states exist, i.e. dim(7i^) = dim(7iB) = 3, is 

H = 2(|00)(00| + |11>(11| + |22)(22|) 

+ |02>(02| + |10)(10| + |21)(21| - 3|^+>(V^+| , (A2) 

where = "^X]i=oN*)- '^^'^ ground-state energy of 
this Hamiltonian is —1, the minimum separable energy 
is and there are bound entangled states with energy 
as low as (3 - 2v^)/3 ~ -0.1547. Although E^^^ = 
the semidefinite program IjAip would return —0.1547 for 
this Hamiltonian, the energy of the minimum energy pos- 
itive partial transpose state. Implementing higher order 
programs as per 45] would give more and more accu- 
rate estimates of the true minimum separable energy, 
-Esop = 0. Furthermore there is a small range of tem- 
peratures, 1.256 < fcsT < 1.271, over which the thermal 
state has energy less than zero, so it is certainly entan- 
gled, but has positive partial transpose. Over this range 
of temperatures the Hamiltonian witnesses the bound en- 
tanglement of the thermal state. 

Examples of Hamiltonians where all low-energy states 
are bound entangled may be constructed from unextend- 



able product bases 48]: a set of product states for which 
the orthogonal complement contains no product states. 
To construct the Hamiltonian we let the unextendable 
product basis span the excited-state manifold, and its or- 
thogonal complement the ground-state manifold. In this 
extreme example, all thermal states with energy within 
the entanglement gap are bound entangled. 



APPENDIX B: ENTANGLEMENT-GAP 
TEMPERATURE OF BIPARTITE SYSTEMS 

In this Appendix, we investigate the entanglement-gap 
temperature of bipartite Hamiltonians. For this purpose, 
we define a completely entangled subspace of a multi- 
partite Hilbert space as one that contains no separable 
states. The antisymmetric subspace of two systems is an 
example. One might wonder whether it is possible to find 
a Hamiltonian with a completely entangled ground-state 
manifold that is larger than the antisymmetric subspace 
so as to achieve a higher entanglement-gap temperature 
than the symmetric projector (|18|l . In ,62j the maximum 
dimension of a completely entangled subspace of many 
parties was investigated: for two c?— dimensional systems 
a basis was given for a completely entangled subspace of 
maximum possible dimension d^ — 2d+l. This subspace 
contains the antisymmetric subspace. 

A natural candidate for a Hamiltonian with a high 
entanglement-gap temperature is thus the Hamiltonian 
with such a subspace at energy zero and its orthogo- 
nal complement at the highest energy. To find its en- 
tanglement gap we could, in principle, use a sequence 
of semidefinite programs as described in Appendix IXI 
However, as the dimension increases we need to imple- 
ment increasingly higher order tests to ensure conver- 
gence and computer memory requirements become pro- 
hibitive. These programs always return a lower bound on 
the entanglement gap. Alternatively, we can bound the 
gap from above by choosing random pure product states 
^ and evaluating their energies. The lowest energy of a 
large number of trial states provides an upper bound on 
the entanglement gap and thus on the entanglement-gap 
temperature. 

Figure |31 compares the behavior of the entanglement 
gap temperature as a function of d for the three Hamil- 
tonians considered above, Hmc = I — \4>d){4'd\, Hs = Us 
and -ffcGs = I — IIccsj where Hccs is the projector onto 
the completely entangled subspace of maximum dimen- 
sion. We see that the entanglement-gap temperature 
of Hccs ~ I — IIccs is generally comparable to that of 
Hnic = I — \<l>d){<i>d\- This result is due to the fact that 



^ To create random pure quantum states we draw the components 
from A/'(0, 1), i.e. a normal distribution with mean zero and vari- 
ance one, and normahze the state. This sampling is equivalent to 
choosing states according to the Haar measure (see for example 
Appendices A and B of Ref . l63l ) . 
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the entanglement gap for -ffccs is quite small, thus re- 
sulting in low entanglement gap temperature despite the 
large ground-state degeneracy. 




FIG. 3: Comparison of entanglement-gap temperature as a 
function of dimension of the subsystems for the three bipar- 
tite Hamiltonians: crosses correspond to -ffmo = I — \(t>d){(l}d\, 
pluses to Hs ~ lis and bounding bars to Hccs ~ I — flcos- 

Another method for constructing completely entangled 
subspaces is as the orthogonal complement of unextend- 
ahle product bases 48] . We have constructed Hamilto- 
nians with completely entangled ground-state manifolds 
from a number of known unextendable product bases 
and have always found entanglement-gap temperatures 
significantly lower than that of the symmetric-projector 
Hamiltonian. 

We thus have good evidence that the symmetric- 
projector Hamiltonian has the highest entanglement-gap 
temperature of Hamiltonians with all energy eigenvalues 
either zero or one. The completely general case where 
there can be intermediate energies as well is beyond the 
scope of this work. 



APPENDIX C: MAXIMUM 
ENTANGLEMENT-GAP TEMPERATURE FOR 
TWO QUBITS 

In this Appendix we investigate the entanglement tem- 
perature of two qubit systems, and provide evidence that 
the Heisenberg antiferromagnetic Hamiltonian has the 
highest scaled entanglement-gap temperature. We scale 
all two-qubit Hamiltonians so that the ground-state en- 
ergy is zero, the maximum energy is one, and there are 
two intermediate energies, < Ei < E2 < 1. The 
antiferromagnet has the singlet at energy zero and all 
triplet states at energy one; its scaled entanglement gap 
and entanglement-gap temperature are gE = 1/2 and 
t£ = l/log,(3). 

We present two lemmas leading to a theorem that any 
Hamiltonian with a maximally entangled ground state 



has an entanglement-gap temperature lower than that of 
the Heisenberg antiferromagnet. 

Lemma 5. Let H and H' be two multipartite Hamilto- 
nians with entanglement gap temperatures Te and T'^, 
respectively. If there is a separable state, pscp such that 



tr[i?'pscp] < U'{Te), 



(CI) 



where U'{T) is the thermal energy of H', then < Tj 



Proof: tr[i/'psop] 
minp,„pe5 tr[H' p^ep] 



is an upper bound on 
By definition U'{T^) 



Sep 

Kep^ SO the 



result follows from the fact that U' (T) is a monotonically 
increasing function of T. □ 

Lemma 6. Any Hamiltonian, H' , with Ei < 1/4 has 
an entanglement- gap temperature less than that of the 
Heisenberg antiferromagnet. 

Proof: We use the fact that, for two qubits, all two- 
dimensional subspaces contains a separable state [6^ . 
Thus, there must be a separable state in the subspace 
spanned by \Eo) and and this separable state must 
have energy less than or equal to Ei. 

We now apply LemmaElwith this separable state pscp- 
Because Ei is the lower of the two intermediate ener- 
gies, the Hamiltonian, iJ", with the same eigenstates and 
eigenenergies as H' , except that E2 = Ei will certainly 
have a lower thermal energy at any particular tempera- 
ture than H', U"{T) < U'{T), V T. The thermal energy 
U"(T) is easily calculated; with it, we find a value of Ei 
that satisfies the condition 

Ei<U"iTE = l/log,{3)) Si < 1/4. (C2) 

Thus, if El < 1/4 then tr[H'pscp] < U"{Te = 
l/log^(3)) < U'{Te = l/log^(3)), so H' has a lower 
entanglement-gap temperature than the Heisenberg an- 
tiferromagnet, as required. □ 

Theorem 4. Any Hamiltonian, H' , with a maximally 
entangled ground state has an entanglement- gap temper- 
ature less than that of the Heisenberg antiferromagnet. 

Proof: Given a Hamiltonian with a maximally entan- 
gled ground state we can use local unitaries to trans- 
form to a Hamiltonian with the singlet as its ground 
state |£:o) = (|0)|1) - |1)|0))/V^. By Lemma □ and 
the invariance of the spectrum under any unitary, this 
Hamiltonian has the same entanglement-gap tempera- 
ture. The excited eigenstates for this Hamiltonian all 
lie in the symmetric (triplet) subspace. We express 
the excited states in their Schmidt decompositions as 
m = A,|0.)|0.) + v/r^|l,)|l,), where i = 1,2,3. 

We present two separable states, one of which has en- 
ergy less than the threshold for any Hamiltonian. The 
first is psep = \ A){A\ ® \B){B\ where 

\A) = (|0i) + |1i))/a/2 , \B) = (|0i) - |li))/%/2. (C3) 
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The energy of this state is at most tr[ff'pscp] = (-£'1 + 
l)/4. For a given Ei this energy will be less than U'{Te) 
for E2 greater than a certain lower bound, E2. i?2^(i?i) 
is defined implicitly by 



tr[i?'p, 



scpj 



U'{TE^l/logM)- 



(C4) 



This equation is transcendental and so it is not possible 
to find an explicit functional form for E!^{Ei). 

The second low-energy separable state that we consider 
ispscp = \A){A\®\B){B\ where \A) = IO3) and |B) = jig). 
The energy of this state is at most tr[H' pscp] = E2/2. 
This energy will be less than U'{Te = l/loge(3)) for E2 
less than a maximum value, Ef{Ei), defined by 



E2/2 = U\TE = l/\og,{3)). 



(C5) 



We numerically solve the two equations, and l|('5p . 
for Ef{Ei) and E2'^{Ei), respectively. From Lemma 
it is only possible that T'j^ > Te = l/loge(3) if Ei > 
1/4. However, it is straightforward to calculate nu- 
merically Ef{Ei) < Ef{Ei) in this region, so that 
for any {Ei,E2), there is a separable state with energy 
less than U'{Te = l/logg(3)). Lemma requires that 
T'e < Te, so the Heisenberg antiferromagnet has the 
highest entanglement-gap temperature of any bipartite 
Hamiltonian with a maximally entangled ground state. 
□ 

A generic two-qubit Hamiltonian has a non-maximally 
entangled ground state, so it is still possible that such 
a Hamiltonian possesses a higher entanglement-gap tem- 
perature than the Heisenberg antiferromagnet. To pro- 
vide numerical evidence that such a Hamiltonian does 
not exist, we generated random Hamiltonians by drawing 
the two intermediate energy levels from a uniform distri- 
bution. Because no bound entangled states exist for two 
qubits, the semidefinite program of Appendix fXlproduccs 
the entanglement gap. We then calculated the entangle- 
ment gap temperature numerically. We generated 10^ 
random Hamiltonians and calculated their entanglement- 
gap temperature in this way. None were found to have 
an entanglement-gap temperature higher than that of the 
Heisenberg antiferromagnet, providing strong evidence 



that it possesses the highest possible entanglement-gap 
temperature. 



APPENDIX D: TRANSVERSE FIELD XY 
MODEL 

The transverse field XY model is defined by the cou- 
pling Hamiltonian Eq. (|34|l . To find the minimum-energy 
separable state |A)|_B) we parameterize the two factors 
as 

|j) =cos0,|T)+e^^^sin^,|i), 3=A,B, (Dl) 

where Q < 9j < 7r/2, < 0j < 2ti. We then calculate the 
energy of the product state |j4)|i?) as a function of the 
four parameters: 

{A\{B\H^]^\A)\B) = |(cos20^ + cos20b) 
+ {^-j^) cos(j)A sin 20^1 COS0B sin20B 

(i^) sin (/)^ sin 26*^ sin 0s sin 26*5 , (D2) 

and optimize over this space to find the lowest energy 
separable state. The result is: 



^=^1 + ^' (D3) 
III), A>l + 7 



where the ± corresponds to j = A, B, with energy 



Sep 



2(1+7) 



-A, 



, A<l+7 
A > 1 + 7- 



(D4) 



Incidentally, by calculating the spectrum of H^g we can 
identify a curve, A^-l-7^ = 1 on which there is a separable 
state in the degenerate ground-state manifold. The en- 
tanglement gap is therefore zero on this curve, and this 
result remains true for the XY model on an arbitrary 
bipartite lattice (with the appropriate magnetic field in 
the coupling Hamiltonian). 
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